Molecular docking analysis of imeglimin and its derivatives with estrogen receptor-alpha

Estrogen receptor-α (ER- α) is a principal endocrine regulatory protein in breast cancer. The progression of ER-α positive breast cancer is slowed by selective estrogen receptor modulators such as Tamoxifen. But, long term therapy with Tamoxifen leads to resistance. Therefore, it is of interest to document the Molecular docking and pharmacokinetic analysis of imeglimin derivatives with ER-alpha. Among the 166 derivatives of Imeglimin, only five derivatives were shortlisted after toxicity testing. The selected derivatives showed good binding affinity with favorable pharmacokinetic profiles. The selected compounds of Imeglimin were found to possess excellent anticancer potential and could be considered as novel, cost-effective anticancer agents effective against ER positive breast cancer for further investigation.


Abstract:
Estrogen receptor-α (ER-α) is a principal endocrine regulatory protein in breast cancer.The progression of ER-α positive breast cancer is slowed by selective estrogen receptor modulators such as Tamoxifen.But, long term therapy with Tamoxifen leads to resistance [1].Therefore, it is of interest to document the Molecular docking and pharmacokinetic analysis of imeglimin derivatives with ER-alpha.Among the 166 derivatives of Imeglimin, only five derivatives were shortlisted after toxicity testing.The selected derivatives showed good binding affinity with favorable pharmacokinetic profiles.The selected compounds of Imeglimin were found to possess excellent anticancer potential and could be considered as novel, cost-effective anticancer agents effective against ER positive breast cancer for further investigation.

Background:
Breast cancer is a significant worldwide health concern.One in every 8 females is predicted to be diagnosed with invasive breast cancer in their lifetime [2].Breast cancer has been treated with a variety of chemotherapeutic drugs, hormonal agents, and targeted drugs.The patient is severely incapacitated by the adverse effect profile of the various drugs classes [3].Furthermore, a significant proportion of the population from lower socioeconomic groups has poor compliance with therapy due to its high cost.The growth and spread of breast cancer, especially hormone receptor-positive breast cancer, are significantly influenced by estrogen receptors (ERs) [4].Apart from the effects on cancer, binding of chemical substances to the estrogen receptor (ER) can lead on to endocrine disruption [5].ER-alpha (ERα) and ER-beta (ERβ) are the two primary forms of estrogen receptors.ERα is the most prevalent form found in breast tissue and is more frequently linked to breast cancer [6].Estrogen receptor-α (ERα) is a 66 kDa protein having a ligand binding domain of 245 residues [7].ER-α is one of the 48 nuclear intracellular receptors [8].It is one of the two major receptors for the endogenous estrogen, 17β-estradiol (E2) [7].ERs are also present in the plasma membrane and mitochondria [9].The transcriptional activities of ER-α can be stimulated by second messenger pathways [10].The altered expression of ERα can be detected in breast cancer using aptamers [11].The progression of ER-α positive breast cancer is slowed by selective estrogen receptor modulators such as Tamoxifen.But, long term therapy with Tamoxifen leads to resistance [1].Moreover, Tamoxifen leads to increased risk of endometrial cancer, stroke and pulmonary embolism owing to its estrogen agonistic action at certain organs [12].Hence, finding new and better drugs is crucial because it is still challenging to address such problems.Recently, Imeglimin, a promising medication for diabetes, structurally similar to Metformin, has drawn attention for its complex pharmacological profile, which includes antiinflammatory and anti-proliferative properties [13].The investigation of Imeglimin and its derivatives in breast cancer is a paradigm-shifting step in pharmaceutical research that could lead to the introduction of a new class of medications for the treatment of this common malignancy.Therefore, it is of interest to document the Molecular docking and pharmacokinetic analysis of Imeglimin derivatives with ER-alpha.

Materials and methods: Ligand preparation:
The structure of the ligand (Imeglimin) was obtained from PUBCHEM database (PubChem CID: 24812808).The chemical structures of various derivatives of Imeglimin were drawn using Chemsketch software.The designed structures were then optimised using Openbabel software.A total of 166 derivatives were designed and the bond lengths and angles were standardized using the clean structure command.In addition, IUPAC name was added using the software [14].

Toxicity Prediction and Pharmacokinetic analysis:
Once the derivatives were developed, simplified molecular input line entry system (SMILES) were created using swissADME, an online tool from Swiss Institute of Bioinformatics (SIB).(http://www.swissadme.ch/)[15].They were screened for various toxicity parameters such as AMES toxicity, Acute oral Toxicity, Carcinogenicity and Rat acute toxicity LD50 using pkCSM.(http://biosig.unimelb.edu.au/pkcsm/)[16].Only those derivatives which exhibited least toxicity were included for further research.Using the pharmacokinetic parameters, the compounds were then scanned for "Lipinski's rule of 5".The shortlisted compounds that did not violate "Lipinski's rule of 5" were then chosen for further study and analysis [17].

Target selection:
The target protein structure was retrieved using RCSB protein databank.Estrogen receptor-alpha (ER-α) (PDB ID: 1A52) was chosen as our receptor target.

Molecular docking:
Protein-ligand interaction analysis is a vital step in virtual screening of potential ligands [18].Discovery studio software was used for initial docking analysis.Receptor ligand interactions were checked and the amino acids that showed significant interaction using hydrogen bond, vanderwaals, etc. were noted (Figure 1).The water molecules and ligands bound to the target molecule (ER-α) were removed and a clean structure was saved as a PDB file.
PyRx software was then used for auto docking.Both the ligand and target molecule were converted into pdbqt format.Using autodock wizard, docking of ligand with the target was performed and the output was obtained.Among the ligand derivatives, the ones with RMSD <3 Å and higher binding affinity were selected.According to the study done by Gonzalez TL et al, RMSDs of ligands docked into human, mouse and rat ER-α were 0.49 Å (human-mouse), 1.19 Å (human-rat) and 0.18 Å (mouse-rat) [19].Using Ligplot+ software, the final docking pose of each derivative with the target protein was obtained (Figure 2).

Results and Discussion:
A total of 166 derivatives of Imeglimin were developed.After toxicity testing, only five derivatives were shortlisted after excluding the toxic derivatives.Figure 3 depicts the 2D structures of the finalized structures of Imeglimin and its derivatives.All the five derivatives of Imeglimin with the Estrogen receptor-alpha (ER-α) showed the binding energies ranging between -6.4 to -7 Kcal/mol which were comparable with that of Tamoxifen (-6.68 Kcal/mol) (Table 1 Acute toxicity studies were performed to analyse the toxicological profile of Imeglimin and its derivatives.The Ames test is a renowned in vitro bacterial mutagenicity assay used to determine the genotoxic potential of substances [29].The negative results obtained from the Ames test for Imeglimin and its derivatives indicate that they possess no mutagenic activity.This implies a minimal risk of inducing mutations in bacterial DNA under the tested conditions.This reinforces the idea that Imeglimin is not genotoxic, which is an important consideration when assessing the safety of intended for therapeutic use.Based on the results of the acute oral toxicity, Imeglimin and its derivatives have an estimated fatal dose for 50% of the tested population (LD50) ranging from 0.4043 to 0.5416.This relatively low LD50 value supports the placement of the compounds into Toxicity Class III.Substances of the toxicity class III are defined as having moderate toxicity, with LD50 values usually falling between 0.1 and 1.0 g/kg body weight.Imeglimin and its derivatives fall into this range, suggesting a moderate degree of acute toxicity.Furthermore, it is reassuring to find out that preclinical research has not shown any carcinogenic characteristics.Carcinogenicity studies provide critical insights into the potential long-term risks associated with exposure to a substance.The lack of carcinogenic effects observed indicates that there is no increased risk of cancer development, further supporting their safety profile (Table 3).
The efficacy and safety of imeglimin derivatives as possible therapeutic agents are largely dependent on their pharmacokinetic properties.The data obtained about their pharmacological characteristics is significant and includes information on their absorption, permeability, plasma protein binding, and inhibition of cytochrome P450 (CYP) enzymes, skin penetration and clearance.Imeglimin derivatives showed high gastrointestinal (GI) absorption which suggests that these compounds are efficiently absorbed from the GI tract into the systemic circulation.This is a favorable characteristic for oral medications, as it indicates a high bioavailability and potential for effective therapeutic action.The absence of blood-brain barrier (BBB) penetration is an important finding because it implies that derivatives of Imeglimin are unlikely to cross the BBB to have an impact on the central nervous system (CNS).This characteristic may lower the possibility of CNS-related side effects, which could improve their safety profile.Plasma protein binding can impact the distribution, metabolism, and elimination of drugs in the body.The observed range of plasma protein binding for Imeglimin derivatives is 7.571% to 22.203%.This is below the optimum plasma protein binding limit of 90%.
High plasma protein binding correlates with low therapeutic index.The lack of inhibition of CYP 1A2, CYP 2C19, CYP 2C9, CYP 2D6 and CYP 3A4 by Imeglimin derivatives is a positive finding, as it indicates that these compounds are unlikely to interfere with the metabolism of other drugs that are substrates for these CYP enzymes.This suggests a low potential for drugdrug interactions involving these pathways.Imeglimin derivatives were found to exhibit a range of skin permeation (-7.23 to -8.28), indicating variable degrees of permeability through the skin.Skin permeation is relevant for topical formulations and can influence the effectiveness of transdermal delivery systems.The range of clearance observed for Imeglimin derivatives (6.015 to 7.513 ml/min/kg) indicates moderate range of clearance from the body.Clearance is an important pharmacokinetic parameter that influences the dosing regimen and overall exposure of an individual to a drug.In summary, the pharmacokinetic parameters of Imeglimin derivatives, including their absorption, permeation, plasma protein binding, inhibition of CYP enzymes, skin permeation, and clearance, provide significant insights into their potential as therapeutic agents.These findings suggest that Imeglimin derivatives have favorable pharmacokinetic profiles, with efficient GI absorption, minimal BBB permeation, low plasma protein binding, no inhibition of major CYP enzymes, and variable skin permeation and clearance rates (Table 4).
Breast .While some research has been done using hepatocellular carcinoma cells, there are no prior studies that have been published in the literature to support Imeglimin's anticancer properties with respect to breast cancer.Hence this molecular docking study was undertaken to ascertain the anticancer effects of Imeglimin and its derivatives in breast cancer.After testing for toxicity, five non-toxic candidates were selected among 166 derivatives of Imeglimin, thus excluding the potentially toxic derivatives.Molecular docking showed that all the five derivatives of Imeglimin formed three hydrogen bonds with the important residues of the estrogen receptor-alpha (ERα).They also exhibited favorable binding energies with ER-α, comparable to that of tamoxifen.The selected derivatives showed a bioavailability score of 0.55, thus exhibiting excellent pharmacokinetic characteristics and substantial systemic absorption.Observation of no violations of Lipinski's rule of five suggests good drug-likeness.
Toxicity studies such as AMES test, acute oral toxicity (LD50) and carcinogenicity were performed.The Ames test revealed no mutagenicity.Based on the LD50 values, Imeglimin and its derivatives were classified into toxicity Class III (moderate toxicity).Preclinical studies proved the absence of any carcinogenic characteristics.Hence these studies provide valuable inputs regarding the safety profile of Imeglimin and its derivatives.Analysis of various pharmacokinetic characteristics showed good gastrointestinal absorption, indicating a high degree of bioavailability, no permeation of blood-brain barrier (BBB), limited plasma protein binding, no inhibition of major CYP enzymes (1A2, 2C19, 2C9, 2D6, 3A4) and moderate range of clearance.In silico techniques are very crucial to establish the 3 R concept-Reduction, Replacement and Refinement.This is useful as an alternative to animal experiments and also reduces the large expenditure involved in research [37].Molecular docking has been a successful approach for determining realistic inhibition mechanisms and ligand-protein interactions.More negative the binding energy (BE), more effective is the ligand binding to the target protein.Spiriti et al have developed a flexible docking approach based on mixed-resolution Monte Carlo (MRMC), to offer a balance with speed, protein flexibility and sampling power [38].Using molecular docking, the selected Imeglimin derivatives exhibited substantial binding (lower BE) with the ER-α comparable with Tamoxifen.This denotes that these compounds may find application as ER-α inhibitors against breast cancer.

Conclusion:
Five non-toxic derivatives of Imeglimin with favorable pharmacokinetic profiles are identified for further consideration as anticancer agents.

Figure 1 :
Figure 1: Amino acids of the protein target Estrogen receptoralpha (ER-α) showing significant bonds and interactions with Tamoxifen generated using Discovery studio software.

Table 2 )
[21]ing affinity of 0.03% as compared to E2[21].All the compounds showed three hydrogen bond interactions each with Asn359, His356 and Arg363 of ER-α with a distance of 2.52, 2.77 and 3.28 Å respectively.Elshal et al. showed that the amino acids Leu A308, Thr C334, Val A368, Thr A371 in the ERα protein interacted with a lectin protein, Concanavalin-Aa.There was a good interaction between Concanavalin-Aa and ERα proving its antagonistic effect on ERα and also showed synergistic action with Tamoxifen in breast cancer [22].In the study by Masand et and maintains the stability of the complex [25].A study done by Patidar et al. tested 40 inhibitors of mTOR receptor protein against breast cancer, among which SF1126 showed the best docking score of -8.705 [26].According to Marwa F Ahmed et al., a series of quinazoline derivatives were synthesized and docking analysis was done against estrogen receptor alpha.All the .According to the study done by Warude et al., indole based benzamides showed only one violation as per Lipinski's rule of 5 and a bioavailability score of 0.55 was obtained [28].